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Abstract 

Charlier configurations provide a combinatorial model for Charlier polynomials. 
We use this model to give a combinatorial proof of a multilinear generating function 
for Charlier polynomials. As special cases of the multilinear generating function, 
we obtain the bilinear generating function for Charlier polynomials and formulas 
for derangements. 

1 Introduction 

Charlier polynomials have been studied using combinatorial methods in [5], [8], [lO], 
^15j . and |il6j|. In this paper, we prove a multilinear generating function for Charlier poly- 
nomials using the combinatorial model of Charlier configurations p^fTT] and the approach 
of Foata and Garsia [6] in their proof of Slepian's multilinear extension of the Mehler for- 
mula for Hermite polynomials [1^. We then obtain some formulas for derangements as 
special cases of this generating function. 

The Charlier polynomials are usually defined by the formula 



c„(a, r) = 2Fo{-n, -a; -r ^) =^i^] 



where {u)k = u{u + 1) ■ ■ ■ {u + k — 1). In order to assign convenient weights in the 
combinatorial model, we work with renormalized Charlier polynomials C„(a, r) defined 
by 



k=0 ^ ^ 
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Our main result is the multilinear generating function 
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where each sum runs over all k x k symmetric matrices {fiij) with non-negative integral 

k 

entries and with diagonal entries zero, = ^ riij for 1 < i < k, and 

To give a combinatorial proof of ([1]), we begin with a discussion of Charlier configu- 
rations and their representation by digraphs in section [2l Then in section |3] we give the 
combinatorial proof of the multilinear generating function. The main idea of the proof is 
to show that both sides of the formula count the same set of digraphs. We discuss the 
special cases of the multilinear generating function in section HI 



2 Charlier Configurations 

Let [n] denote the set {1, 2, . . . , ra}. 

Definition 2.1 A Charlier configuration on the set S is a pair $ = ((A, cr), i?), where 
{A^B) is an ordered partition of S and a is a permutation of A. 

A Charlier configuration is called a partial permutation in [10]. The configuration $ 
can be represented by a digraph with vertex set S and with an edge from i to j if and 
only if a{i) = j. Figure [1] shows a Charlier configuration on [10]. 
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Figure 1: Charlier configuration on [10] 



Here A = {1, 3, 4, 5, 6, 7, 9, 10}, B 
cycle notation. 



{2,8}, and a = (7) (4 9 6) (1 5 10 3) in disjoint 



2.1 Combinatorial interpretation of Charlier polynomials. 

We assign a weight to a Charlier configuration $ by assigning a weight a to each cycle 
of o and a weight r to each point of B. Then the weight of $ is a'^^'^'^'^V^' where cyc(o") 
denotes the number of cycles in a. (If we had not renormalized the Charlier polynomials, 
we would assign the weight —a to each cycle of a and the weight 1/r to each point of 
A.) We use the following well-known facts about generating functions for permutations, 
which are proved, for example, in [131 p. 19]. 

Fact 1. Tl=oc{n,k)a}' = (a)„ where c(n, k) is the number of permutations of [n] with 
exactly k cycles (the unsigned Stirling number of the first kind). 

Fact 2. The exponential generating function for all permutations, with cycles weighted 
by a, is (1 — z)~"'. 

Let ^5 denote the set of Charlier configurations on S. Then it follows easily from 
Fact 1 that Cn{a, r) is sum of the weights of the elements of ^[n\- 



3 Combinatorial Proof of the Multilinear Formula 

We assume that the reader is familiar with enumerative applications of exponential gener- 
ating functions, as described, for example, in Chapter 5] and [3J. The product formula 
and the exponential formula for exponential generating functions discussed in these ref- 
erences play an important role in the combinatorial proof of the multilinear formula. The 
theory of species (as used in (TU] and [TT]) could be used to provide a proof of the formula 
as well. 

The formula ([1]) could be proved by interpreting it as a multivariable exponential 
generating function in the variables Xjj, which would require the use of digraphs with 
multiple sets of labels. The proof is simpler if we use exponential generating functions in 
only one variable, so that we can use a single set of labels. To accomplish this, we rewrite 
the formula by replacing Xij with zxij . Now we can think of the formula as an exponential 
generating function in the single variable z. The formula is now 
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We will prove this formula, which is equivalent to ([T]). We begin with a description of the 
digraphs counted by the left side of the formula. 
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3.1 Digraphs counted by the left side. 

We can rewrite the left side of ([2]) as follows: 

n>0 ■ (n,,) ^J} l<i<j<k 
Z, nij=n 

i<j 

Let {Nij)i<i^j<k be an ordered partition of [n] such that \Nij\ = riij. For j > i, let 
Nji = Nij. Let iVj = Uj^iNij. Then A^j fl A'j- = Nij. Since rij = Ylj'^^ij^ follows that 
|A^j| = rii. Let if be the set of all ordered tuples {{Nij), <I>i, . . . , such that 

1. {Nij) is an ordered partition of [n] with the above properties. 

2. Each $i is a Charlier configuration on Ni, i.e., $i G ^at.. 

Then each point of [n] is in exactly two configurations. This follows from the fact that 
each point is in exactly one Nij and Ni fl Nj = Nij. To the Charlier configuration 
$j = {{Ai, ai), Bi) we assign the weight a^^'^'^'^Vj'^'. We also assign an additional weight of 
Xij to each point of Nij. The weight of a tuple in H is defined to be the product of the 
weights of its constituent Charlier configurations and its points. Then it is easy to see 
that the left side of ([2]) is the exponential generating function for H with these weights. 

We associate a digraph to a tuple {{Nij), $i, . . . , in H by superimposing the di- 
graphs of these k Charlier configurations on [n] in which each <l>j is on of these vertices. 
Figure [2] shows such a digraph for k = 3. The configurations <l>i,$2,'^'3 are respectively 
represented by sohd lines, dashed lines, and dotted lines. Each vertex is in exactly two 
configurations and this is indicated by the two different circles around each vertex. 

The tuple and the configurations corresponding to Figure [2] are given by 

Ni2 = {5, 8, 9, 11, 14, 16, 17}, Ari3 = {1,2,6,12,15}, and iVas = {3,4,7,10,13,18}. 
$1 = {{Ai,ai),Bi), where 

Ai = {1, 8, 9, 11, 12, 15, 16}, Bi = {2, 5, 6, 14, 17}, 

(Ti = (9 16 15) (1 12 11 8). 
$2 = ((^2, cr2), -B2), where 

A2 = {3, 4, 5, 7, 8, 10, 11, 17, 18}, B2 = {9, 13, 14, 16}, 

a2 = (5 17 8 11) (3 10 7 4 18). 
$3 = {{A3,a3,),B3), where 

As = {1, 4, 6, 7, 10, 13, 18}, ^3 = {2, 3, 12, 15}, 

as = (6 13) (1 10 7 18 4). 

We may identify H with the set of these digraphs, for which so that the left side of ([2]) 
is the exponential generating function. We now enumerate these digraphs in another way: 
We consider their connected components, which are of three types, and use the product 
formula for exponential generating functions to show that the right side of ([2]) is also a 
generating function for H. 
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Figure 2: Digraph counted by left side for /c = 3 



3.2 Connected components of digraphs in H. 

Let {{Nij), $1, . . . , ^k) be a tuple in H, where = ((A, en), Bi) G C„,, for 1 < i < k. 
The connected components of the digraph representing this tuple are of the following 
three types: 

For i < j, a type ly connected component is an isolated vertex which is in $j and 
$j but not in or (Xj. In Figure [21 vertex 2 is of type I13 and vertex 14 is of type I12. 
Such a vertex belongs to BidBj and is weighted by riTjXij. It follows that the exponential 
generating function for digraphs all of whose components are of type Ijj, which we call 
type lij digraphs, is e^i^j^ij^, 

A type 2j connected component is a cycle of ai in which no vertex is in any 
other aj. In Figure [2], the cycle (6 13) is a type 2^ component and (9 16 15) is a type 2i 
component. The cycle of a type 2j component weighted by and each vertex of the 
cycle is in some Bj and so is weighted by and Xij . A type 2j digraph is a digraph 
in which every connected component is of type 2j. Such digraphs can be considered as 
permutations in which each cycle is weighted by and each vertex is weighted by some 
Tj and Xij for some j. By a slight modification of Fact 2 in subsection 12. 1^ it follows that 
the exponential generating function for type 2j digraphs is (1 — 2; ^ rjXij)~°'\ 

Any connected component that is not of type Ij or 2ij is called a type 3 connected 
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component. In a type 3 connected component, every vertex is in at least one permutation 
and every cycle contains at least one vertex that is also in another permutation. 

The type 3 connected component from Figure [2] is shown in Figure [31 A type 3 




Figure 3: Type 3 connected component 



digraph is a digraph all of whose connected components are of type 3. We say that a 
type 3 digraph is reduced if every vertex is in two permutations. Thus, a reduced type 3 
digraph on n vertices is an ordered tuple {{Nij), $i, . . . , where each $j = ((A, o'i), 0); 
i.e., each $j is simply a permutation ai on Ui vertices. Figure H] shows a reduced type 3 
digraph on 7 vertices. Since each cycle of cTj is weighted by a,, the exponential generating 




Figure 4: Reduced digraph on 7 vertices 
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function for reduced type 3 digraphs is 



n>0 (mA -.^-V-^, ^ l<i<k l<i<j<k 
— \ ijj l<i<]<k — J — 

2. nij=n 

i<j 

Each vertex in a reduced type 3 digraph has two outgoing edges belonging to two different 
permutations. Any type 3 digraph can be obtained from a reduced type 3 digraph by 
replacing each outgoing edge at every vertex by a sequence of ordered edges. An outgoing 
edge in cTj is replaced by a sequence of edges, such that each new vertex is in cXi, but not 
in any other aj. Hence each new vertex is weighted by rj and Xij for some j. Thus the 
exponential generating function for type 3 digraphs is 



7 . rt 1 

X ^ ' t' . ■ ■ \ ^'1 I ri ^ mm • 

X ■ 



It follows from the product formula for exponential generating functions that the ex- 
ponential generating function for digraphs in H is the product of the generating functions 
for all of the types of digraphs described above and this is 

n nvjx.jz TT \ 
J- J- (l-zy2riX,,Y' 

\ — *' T — r i.^iJ'^i T — r 

^ 2^ fr ^ 11 fi - z"p r-x--Y' "^'^ ' 

n>0 ■ ("ij) l<j<fe ^ f:^. ^ l<i<j<k 



l<i<j<k l<i<k ... 
-J- 



^ l<i<7'<fc jj^i 



which is equal to the right side of ([2]) 



4 Specializations 

For k = 2, the only parameter riij in ([1]) is riu, and rii = n2 = ni2. If we write n for ^12, 
a for ai, 6 for 02, r for ri, s for r2, and a; for X12 then the multilinear formula ([1]) reduces 
to the bilinear formula 

j:C,(aM.fl ^ e--^ _ (3) 

n>0 n>0 ^ ' ^ ' 

Formula (2.47) in Askey's book ^ is equivalent to the case of ([3]) in which a and h are 
negative integers, and the general case is easily derived from this. Note that a and h are 
switched on the right side of the formula in the book. 
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Similarly, the case = 3 of ([T]) may be written 



I rtiiTL 



^-^ l\ ml n\ 

l,m,n 

{a)l+mib)i+n{c)m+n X^y^Z^ 

^ (l- sx - -rx - tzy+"-+^{l -ry - l\ ml n\' ' ' 

Some special cases of these formulas are worth mentioning. Setting a; = in (jlj) gives 

y^ Cm{a, r)Cn{b, s)C„+„(c, t) — — 
^-^ ml nl 



^ ^{ry+sz)t {0')m{b)n{c)m+n y^Z^ 

~ mn^^~ ty)"^+''{l - tzY+^{l -ry - mini' ^ ^ 

Formula may be viewed as a Charlier polynomial analogue of a formula of Carlitz [1] 
for Hermite polynomials, which is a special cases of Slepian's multilinear extension of the 
Mehler formula [12]. 

By applying the fact that Cn(0, 1) = 1, we can find other simplifications. Thus setting 
6 = and s = 1 in ([3]) gives the usual exponential generating function for Charlier 
polynomials, 

Y^C^(a,r)^ = e^^{l-x)-\ 
^ — ^ nl 

n>0 

Setting c = and t = 1 in (jlj) gives a generalization of ([2]): 

J2 G+„.(a,r)a+„(6,.)|^^ = e-+'^^+- E (i _ _ y^Hl- rx - zy-^^ll' 



Lm,n 



4.1 Permutations. 

The Charlier polynomials can be normalized in another way so as to count permutations 
by cycles and fixed points. We define polynomials Dn{(y,u) by 

where is the set of permutations of [n], cyc~^i(7r) is the number of cycles of vr of length 
greater than 1, and fix(7r) is the number of fixed points of vr. 

We can express the polynomials Dn{<y,u) in terms of Charlier polynomials. To a 
Charlier configuration $ = {{A,a),B) on [n] we may associate the permutation vr of [n] 
such that 71(1) = a{i) for i & A and 7r(i) = i for i E B. Conversely, given a permutation 
71 of [n], the corresponding Charlier configurations may be constructed by choosing an 
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arbitrary subset B of the set of fixed points of vr and taking a to be the restriction of vr 
to A = [n] \ B. This construction yields the relation 

Cn{a,r) = Dn{a,a + r) 

and thus 

Dn{a, u) = Cn{a, u — a). 

So formulas Q-Q may rewritten as generating functions for the polynomials Dn{<y,u). 
Of particular interest are the specializations -Dn(a) = Dn{a,0) = Cn(a,— a), which 
count derangements (permutations without fixed points) by cycles, and D„ = -D„(l) = 
C„(l, — 1), the number of derangements of [n]. 

Making the appropriate substitutions in Q gives 

n=0 n=0 V ' ^' y V ' J 

Formula ([7]) was proved by Gessel [7] as a special case of a generating function for 3 x n 
Latin rectangles, using an approach similar to that of this paper. (See also [Ij and [17J.) 
The corresponding specialization of ([5]) is 

ym^n 



D^{a)Dn{l3)Dm+n{l) 

a)m(/3)n(7)m+n V"" Z 



m,n=0 



and of (ED is 
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{l + (3x- y)'+"(l + ax- zY+l^ n 
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